PROCEEDINGS OF THE 31 st ICRC, LODZ 2009 



1 



Search for High Energetic Neutrinos from Supernova Explosions 

with AMANDA 

Dirk Lennarz*, Jan-Patrick HuIB* and Christopher Wiebusch* for the IceCube Collaboration^ 

*///. Physikalisches Institut, RWTH Aachen University, 52056 Aachen, Germany 
^See http://www.icecube.wisc.edU/collaboration/authorlists/2009/4.html for a full list of authors 



Abstract. Supernova explosions are among the most 
energetic phenomena in the known universe. There 
are suggestions that cosmic rays up to EeV energies 
might be accelerated in the young supernova shell 
on time scales of a few weeks to years, which would 
lead to TeV neutrino radiation. The data taken 
with the AMANDA neutrino telescope in the years 
2000 to 2006 is analysed with a likelihood approach 
in order to search for directional and temporal 
coincidences between neutrino events and optically 
observed extra-galactic supernovae. The supernovae 
were stacked in order to enhance the sensitivity. A 
catalogue of relevant core-collapse supernovae has 
been created. This poster presents the results from 
the analysis. 
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I. Introduction 

Almost a hundred years after their discovery, the 
acceleration mechanisms and sources of the cosmic 
rays remain an unsolved problem of modern astronomy. 
Neutrino astronomy can be an important contribution to 
the solution of this problem. Young supernovae in con- 
nection with a pulsar have been proposed as a possible 
source of cosmic rays with energies up to the ankle. This 
pulsar model can be directly tested by measuring high 
energetic (TeV) neutrino radiation on time scales of a 
few weeks to years after the supernova HU3. 

The AMANDA-II neutrino telescope is located in the 
clear ice at the geographic South Pole and was fully 
operational since 2000. It reconstructs the direction of 
high energetic neutrinos by measuring Cherenkov light 
from secondary muons. The main background are muons 
and neutrinos produced in air showers in the atmosphere. 

This analysis uses 7 years of AMANDA data taken 
during the years 2000-2006 with a total live-time of 1386 
days. The data reconstruction and filtering is described in 
[ 3 1 and the final event sample contains 6595 events. The 
contamination of mis-reconstructed atmospheric muon 
events is less than 5% for a declination greater than 5°. 

II. Pulsar Model 

Rotational energy liberated by a pulsar can be con- 
verted into the energy of relativistic particles [1|. Sec- 
ondary particles, for example pions, are created in the 
interaction with the expanding supernova envelope and 
decay into neutrinos and other particles [1|. In this 
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Fig. 1: Typical supernova model neutrino light curve 

analysis the pulsar model as described in [2] is used. 
Thermonuclear supernovae have no pulsar inside the 
envelope and are therefore not considered by this model. 

The phase of powerful, high energetic neutrino emis- 
sion is limited by two characteristic times: the time 
at which the pion decay time becomes less than the 
time between two nuclear collisions (t w ) and the time at 
which the density of the envelope is sufficiently small 
for accelerated particles to escape into the interstellar 
space without interaction (t c ). The supernova neutrino 
luminosity as a function of time (model light curve) is 
given by: 

■\L (l + ^\ , (1) 

where A is the fraction of the total magnetic dipole 
luminosity Lq (in erg/s) that is transferred to accelerated 
particles and r the characteristic pulsar braking time. 

The shape and length of the model light curve depend 
on the supernova envelope mass (M c ), uniformity (de- 
scribed by a parameter called £) and expansion velocity 
(V), the pulsar braking time and the maximum pion 
energy. An E~ 2 neutrino energy spectrum is assumed 
with an energy cutoff at 10 14 eV. Fig. Q] shows a typical 
model light curve for (,Ri8x 10 3 s and t c « 2 x 10 6 s. 
These values are obtained by choosing M c = 3M Q , 
£ = 1, V = 0.1c and r = lyear. 
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III. Supernova Catalogues 

For this analysis a catalogue of supernovae was cre- 
ated. It combines three different electronically avail- 
able and regularly updated SN catalogues [4|[5||6|. 
A comparison of the three catalogues revealed some 
inconsistencies in the listed information. A consistent 
selection was made with special attention to the objects 
mistaken for a supernova observation, the total number 
of supernovae and the supernova positions. 

Fig. 12 shows the distribution of the 4805 supernovae 
observed between 1885 and 2008. Out of these only 
about 700 supernovae inside the AMANDA data taking 
time are relevant. The clearly visible structure around 
the celestial equator are supernovae found by the Sloan 
Digital Sky Survey-II supernova survey. The nearest and 
best visible supernova for AMANDA was SN2004dj in 
NGC 2403 at a distance of approximately 3.33 Mpc 
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Fig. 2: Distribution of observed supernovae in equatorial 
coordinates with the galactic plane indicated as dashed 
line. Due to the background from atmospheric muons 
only supernovae in the northern hemisphere are relevant. 

This analysis searches for directional and temporal co- 
incidences between neutrinos and supernovae. Therefore 
additional input has to be quantified for each supernova. 
Firstly, the expected neutrino flux has to be determined 
from an accurate distance. The supernova distance can 
be identified with the distance to the host galaxy and can 
be estimated from the redshift. The redshift estimate is 
replaced by a measured distance (e.g. Cepheid variables 
or Tully-Fisher relation) if available. This improves the 
distance accuracy for nearby supernovae, which are most 
relevant. 

Secondly, the explosion date is needed for the temporal 
correlation, but only the date of the optical maximum 
or the discovery date is available. From some well 
observed SNe (e.g. 1999ex and 2008D) it is known that 
the optical maximum occurs around 15-20 days after 
the explosion, which is used as a benchmark. Fig. [3] 
shows the difference between the date of discovery and 
the date of maximum for those cases where the light 
curve was fitted to a template and the date of maximum 
extrapolated backwards in time or found on old photo 



plates. The majority of the supernovae are discovered 
within 20 days after the optical maximum. Hence, the 
discovery is assumed to be typically 20 days after the 
optical maximum. The uncertainty of the explosion date 
is accounted for in the likelihood approach (see next 
section). 
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Fig. 3: Number of days between the optical maximum 
and the discovery if the supernova was discovered after 
the maximum. A linear one day binning is shown on an 
logarithmic x-axis. 

Thirdly, for the individual supernova the needed input 
for the pulsar model is not available. Therefore, all 
supernovae are treated equally. The influence of the 
model light curve on the analysis is tested by defining 
two additional sets of parameters which result in light 
curves with very short and long neutrino emission. 
Hence, altogether three different model light curves 
(typical, short, long) are used. The width of the plateau 
that can be seen in Fig. [T] is 12 days for the typical, 1 
day for the short and 76 days for the long light curve. 
The most realistic assumption for the supernovae in the 
catalogue is that they have individual realisations of the 
parameters of the pulsar model and therefore individual 
light curves between the extreme cases. 

IV. Likelihood Approach 

A new likelihood approach was developed for this 
analysis [7|. Its principal idea is to compare all neutrino 
events from the experimental data sample to every rele- 
vant supernova and evaluate the likelihood ratio (LHR) 
between the hypothesis that this event is signal and 
the hypothesis of being background. This yields a large 
value for a good and small value for a bad match. The 
LHR for all events is summed in order to obtain a 
cumulative estimator, called Q: 

v- Esnp(«I sn M sn ) 

events 



Q 



p(a\BG)p(BG) 



(2) 



where a are characteristic observables of the event. 

The advantage of this likelihood definition is that it 
can be extended to a stacking analysis. Q automatically 
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assigns a small weight to irrelevant combinations of 
neutrinos and supernovae, while relevant ones receive a 
larger weight. Thus, all supernovae from the catalogue 
can be used in the analysis and no optimisation on the 
number of sources is needed. Q is a sum of likelihood 
ratios and therefore its absolute value contains no phys- 
ical information. 

The probabilities in the likelihood sum are constructed 
from properties of AMANDA, the experimental data 
sample and the considered model light curve. p(BG) is 
the probability to have background and is an unknown 
but constant factor. This probability is eliminated by 
redefining Q to Q -p(BG). 

p(a\BG) is the probability that, assuming an event is 
background, it is observed at its specific time and from 
its specific direction. It is factorised into a temporal 
and an angular part. The temporal part corresponds to 
the AMANDA live-time. However, it cancels out with 
the corresponding temporal part of p(a|SN). AMANDA 
does not distinguish between signal and background 
neutrinos and was obviously taking data when the event 
was measured. The angular probability is constructed 
with the normalised zenith angle distribution of the 
experimental data sample (see Fig. 2 in The azimuth 
probability is constant, because AMANDA is completely 
rotated each day and the azimuth is randomised for the 
relevant time scales of this analysis. 

The supernova signal probability consists of p(a|SN) 
and p(SN). The first part depends on the specific event 
and is the probability that an event from a supernova is 
observed at a given time from a given direction. p(SN) 
is the probability to observe a signal from that supernova 
and is estimated for each supernova. 

For p(a|SN) two terms are considered. p(4 f |SN) is 
the probability that a neutrino from a supernova is 
reconstructed with an angular difference ^ relative to 
the supernova direction. This probability is calculated 
from the point-spread function, which is obtained from 
Monte Carlo simulations. The second term p(t, tg-^\SN) 
yields the probability that a neutrino arrives with a time 
offset t — ^sn fr° m me explosion date. This probability 
is taken from a likelihood light curve. 

In order to be less model dependent three generic 
likelihood light curves (typical, short, long) are con- 
structed. They are inspired by the model light curves and 
constructed conservatively in order to not miss signal by 
accidentally looking too early or too late. Hence, if the 
date of the optical maximum is known, the starting time 
for the likelihood light curves (t = 0) is defined to be 
30 days earlier. In case only the date of the discovery 
is known a 50 days earlier starting time is used. This 
makes sure that the explosion is not missed, because the 
time shift to the explosion date is overestimated by about 
15 days for the optical maximum and up to 35 days for 
the date of discovery. The likelihood light curves consist 
of a half Gaussian for t < 0, a plateau for t > and 
another half Gaussian after the plateau. The length of the 
plateau is the full width at 90% of the model light curves 



and enlarged by the uncertainty of the explosion day. 
This uncertainty is bigger if only the date of discovery 
is known. The width of the Gaussian after the plateau 
is the full width at half maximum (FWHM) after the 
plateau of the model light curves. Fig. |4] shows the 
typical likelihood light curve for the date of discovery. 
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Fig. 4: Typical likelihood light curve 

p(SN) depends on the supernova neutrino luminosity, 
distance and direction. The absolute value of p(SN) is 
determined by the supernova neutrino luminosity and is 
a free parameter of this analysis. However, the absolute 
normalisation is not required, because a constant factor 
results in a rescaling of Q and hence only relative values 
are important. All supernovae are assumed to have the 
same neutrino luminosity at source. p(SN) decreases like 
the flux with the square supernova distance. AMANDA 
is not equally sensitive to neutrinos from all directions. 
Therefore the angular acceptance for different supernova 
directions is taken into account. 

V. Signal and Background Simulation 

Q distributions for signal and background simulations 
are used to construct confidence belts with the Feldman- 
Cousins approach to the analysis of small signals (8). 
Each simulated data sample contains 6595 signal or 
background events like the experimental data set. Back- 
ground events are simulated with the zenith angle dis- 
tribution of the experimental data and the AMANDA 
live-time. For the signal simulation a model light curve 
and the AMANDA angular and temporal acceptance is 
simulated. The angular acceptance includes a random 
simulation of assumed systematic uncertainties of the 
measured rate of high energetic muon neutrinos |3|. 
The simulation of the temporal and angular acceptance 
reduces signal events from days with low live-time or 
unfavourable supernova directions. 

The confidence belts are used to estimate the sen- 
sitivity of the analysis. The sensitivity for the long 
model light curve is not compatible compared to Q. 
Furthermore, if the supernovae have short model light 
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curves, the sensitivity is comparable for the short and 
typical likelihood light curves. The typical pulsar model 
is best detected with the typical likelihood light curve. 
Therefore the experimental data is analysed with the 
typical likelihood light curve, because it can cover a 
larger range of possible parameters. 

VI. Experimental Result 

Analysing the experimental AMANDA data with the 
typical likelihood light curve yields: 

Qf y x p P cal = 0.0059 . (3) 

Fig. [5] shows this value in a Q distribution for back- 
ground only. The p value of obtaining a Q value equal 
or bigger than 0.0059 is 73.0%. Hence, the Q value is 
consistent with background and no deviation from the 
background only hypothesis is found. Therefore upper 
limits for the three model light curves are derived. 
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Fig. 5: Q value from analysing the experimental data 
sample with the typical likelihood light curve (horizontal 
line) in a distribution for background only. 

90% upper limits on the signal strength are derived 
from the Feldman-Cousins confidence belts. With the 
help of the signal simulation explained above this value 
converts to the sum of neutrinos from all supernovae. 
Assuming the above model ranking of sources and the 
stacking this limit can also be interpreted as a limit on 
the number of neutrinos from SN2004dj. Tab. U shows 
the obtained upper limits for the typical, short and long 
pulsar model light curve. 



Pulsar model 


All SNe 


SN2004dj 


Typical 


< 5.4 


< 1.0 


Short 


< 4.1 


< 0.9 


Long 


< 67.3 


< 5.9 



The event numbers can be converted to a flux by 
integrating the AMANDA neutrino effective area with 
the expected signal energy spectrum (E~ 2 spectrum with 
cutoff at 10 14 eV). Assuming the typical pulsar model 
for all supernovae and taking the average of the effective 
area over all directions, the 90% upper limit on the flux 
from all supernovae for the plateau of powerful neutrino 
radiation (12 days) is: 

% ■ < 5.2 x 10- 6 ^ . (4) 
a E cm^s 

Using the effective area for the direction of SN2004dj, 
the corresponding 90% upper limit for SN2004dj is: 

-GeV 



dE 



< 8.4 x 10" 



(5) 



TABLE I: 90% upper limits on the number of neutrinos 
from all supernovae and from SN2004dj. 



These limits are valid in the energy range from 1.1 
TeV to 84.0 TeV 

Assuming that the energy range of the pulsar model 
as described in can be extended to higher energies, 
the limits would improve by about 30% and are then 
valid in the energy range from 1.7 TeV to 2 PeV. 

VII. Conclusion 

For the first time the neutrino emission from young 
supernova shells was experimentally investigated. In 
the context of the pulsar model no deviation from the 
background only hypothesis was found. 

For a galactic supernova the expected flux from the 
pulsar model should be sufficient to be detectable by 
IceCube, the AMANDA successor. The sensitivity of 
this analysis might be enhanced by using an energy 
estimator in the likelihood and the individual event re- 
construction error instead of the energy averaged point- 
spread function. 
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